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ABSTRACT 

We provide an exhaustive analysis of the Integrated Sach-Wolfe effect (ISW) in the context of 
coupled Dark Energy cosmologies where a component of massive neutrinos is also present. We focus 
on the effects of both the coupling between Dark Matter and Dark Energy and of the neutrino mass 
on the cross-correlation between galaxy/quasar distributions and ISW effect. We provide a simple 
expression to appropriately rescale the galaxy bias when comparing different cosmologies. Theoretical 
predictions of the cross-correlation function are then compared with observational data. We find 
that, while it is not possible to distinguish among the models at low redshifts, discrepancies between 
coupled models and ACDM increase with z. In spite of this, current data alone seems not able to 
distinguish between coupled models and ACDM. However, we show that upcoming galaxy surveys will 
permit tomographic analysis which allow to better discriminate among the models. We discuss the 
effects on cross-correlation measurements of ignoring galaxy bias evolution, b(z), and magnification 
bias correction and provide fitting formulae for b(z) for the cosmologies considered. We compare three 
different tomographic schemes and investigate how the expected signal to noise ratio, snr, of the ISW- 
LSS cross-correlation changes when increasing the number of tomographic bins. The dependence of 
snr on the area of the survey and the survey shot noise is also discussed. 

Subject headings: (cosmology:) cosmic microwave background; cosmology: miscellaneous; cosmology: 
observations; cosmology: theory; cosmology: large-scale structure of universe 



1. INTRODUCTION 

Several observations made over the recent years, re- 
lated to a large extension to Large Scale Structures 
(LSS) and anisotropies of the Cosmic Microwave Back- 
ground (CMB) as well as the magnitude-redshift relation 
for type la Supernovae have given us a convincing pic- 
ture of the energy and matter density in the Universe 
(iPerlmutter et al 111 9991: iRiess et al 111998b ISt>ergel et al I 
20031 iTegmark et al 1120041 lLarson et al II2O1O0 . 

Baryonic matter accounts for no more than 30% of the 
mass in galaxy clusters while the existence of a large 
clustered component of Dark Matter (DM) seems now 
firmly established, although its nature is still unknown. 
However, they contribute to the total energy density of 
the Universe with only a few percent and about 25% 
respectively. 

No more than another few percent could be accounted 
for by massive neutrinos, but o nly in the most favorable , 
but unli kely case. Acco r ding t o iKristiansen et al I ()2007l ) 
(see also lElgarOv et all (|2002l )) the total mass of neutri- 
nos cannot exceed the li mit of 1.43 eV (see, how ever, 
lLa Vacca efafl (|2009allbl ): IKristiansen etail (12010ft for a 
recent analysis on neutrino mass limits in coupled dark 
energy models). A very small part (10~ 4 ) of the total 
energy density is due to massless neutrinos and CMB 
radiation. 

The model suggested by observations is only viable if 
the remaining 75% is ascribed to the so-called Dark En- 
ergy (DE) responsible for the present day cosmic accel- 
eration. 

Although strongly indicated by the observations, the 
existence of DE is even more puzzling than DM. It can 
be identified with a cosmological constant A or with a yet 



unknown dynamical component with negative pressure. 
On the other hand, its manifestation can be interpreted 
as a geometrical property of the gravity on large scales 
resulting from a failure of Gene r al Re lativity (GR) on 
those scales (see iCopeland et ail (|2006[ ) for a review). 

Within the context of GR, as an alternative to the 
cosmological constant, DE is usually described as a 
scalar field 4>, self-interacting through a suitable po- 



tential V((t>), or a cosmic flu id with negative p ressure 
Peebles fc Ratra I (12003ft / s ee iKoivisto fc MotaT (T2 006) ; 



Koivi sto. Mota fe Pitroul (j2009ft ; iLi, Barrow fc Mota I 
(|2007ft for alternatives and references therein). 

Scalar fields naturally arise in parti cle physics. Fur- 
thermore, if they are tracker fields (jSteinhardt et all 
11999ft . fine tunings associated to the small value of the 
present DE energy density can be significantly alleviated 
unlike the cosmological constant case. 

In addition to self-interaction, a scalar field can in 
principle be coupled to any other field present in na- 
ture. However, in order to drive the cosmic accelera- 
tion, its present time mass is expected to be, at least 
on large scales, ~ H ~ 10~ 33 eU (H being the 
present Hubble parameter). Such a tiny mass gives rise 
to long-range interactions which could be tested with 
fifth-force type experiments. Couplings to ordinary par- 
ticles are strongly constrained by such a kind of exper- 
iments but limits on the DM coupling are looser (con- 
straints on coupling fo r specific models were obt a ined i n 
Maccio' et ail (12004ft : lAmendola fc Quercellinil (2 003) ; 
Olivares et ail (120051): iLee et al I d2006ft : IGuo et al I 
' 2007ft : iMainini fc Bonometto I (|2007ft from CMB, N- 



body simulations and matter power spectrum analysis). 

A possible common ori gin of DM and DE and/or 
a their direct coupling (jWetterich I 119951 : lAmendoIal 



1999; iGasperini et al I 



2000; iChimento et al I 
Farrar fc Peebles I 120041 ). 



2002; Bart olo fe Pietroni I 
12003 IRhodes et al ] [2003; 
would ease one of the 
most critical problems in modern cosmology, the so- 
called coincidence problem: why expansion started to 
accelerate just at the eve of our cosmic epoch, after 
decelerating during all epochs after inflation? Why DE 
and DM have similar densities just now? Because of 
the coupling DM and DE densities keep similar values 
during a long period and the only peculiar feature of the 
present epoch is the recent overtaking of DM density by 
DE density. 

If present, DM-DE coupling could have a relevant role 
in the cosmological evolution affecting not only the over- 
all cosmic expansion but also modifying the DM particles 
dynamics with relevant consequences on the growth of 
the matter density perturbations in both linear and 
nonlinear regime (e.g., on halo de nsity profiles, cluster 
mass f u nction and its e v olution , se e I Wang fe Steinhardt I 
(fl998l): IMainini et al I ti2003aHbf); iKlypin et aiTT2 003): 
iDolag et al I (120041) ;lMaccio' et al | (120041) ;lPerrotta et al I 
IS) 



20041): IMota fe Barrow I (120041): loiivares et all (|2006|) : 



Nunes et al 



20051): IMota fc van de Bruck I 12004) 
.Maor fc Lahav I (120051) : jNunes et al I (120051): iWang 
(120061) :lManera fc Motal (120061): iNunes fc Motal (120061) 



iDutta fc Maor 



. (120071) : IMota et al I (120071) : I Mainini I 
20081 ) : iShaw fc Motal (I2008D : iMotal 120081): IMainini 
2009); IBaldi et all (120091 ): iWintereest fc Pettoriiiol 
( 20101 ): IBaldi fc Pettorino I ()2010l ). LSS is then a 
powerful probe of DE nature which permit to put 
significant constraints on DE parameters. Constraints 
often become even more stringent when data from other 
probes are simultaneously taken into account. 

CMB is another powerful probe of DE nature. In prin- 
ciple, by joining anisotropy and polarization data, DE 
parameters can be significantly constrained. CMB and 
LSS probe the universe at different epochs and are there- 
fore complementary to each other. Future data from high 
resolution CMB experiment like PLANCK Q and LSS 
surveys (EUCLID El, LSST0, DES0, JDEM0, etc) will 
allow to constrain DE up to an unprecedented accuracy. 

In this paper we will focus on the Inte grated Sachs 
Wolfe (ISW) effect (ISachs fe Wolfe I fl967l) . ISW effect 
is a secondary anisotropy of the CMB and a direct sig- 
nature of DE. The effect arises when a photon from the 
last scattering surface passes through a time-dependent 
gravitational potential changing its energy so that addi- 
tional temperature anisotropics are generated. Decay of 
gravitational potentials may occur through cosmic cur- 
vature, in the presence of DE or in alternative gravity 
models. 

Assuming General Relativity is the correct theory of 
gravity and that the Universe is spatially flat, large-scale 
gravitational potentials do not evolve significantly in the 
matter era. Cosmic acceleration, however, causes the 
gravitational potentials to decay making the ISW effect 
highly sensitive to the presence of DE. 

Though difficult to detect directly in the CMB, ISW 



1 http:/ /www.rssd.esa.int/index.php?projcct=planck 

2 http://www.euclid-imaging.net/ 

3 http://www.lsst.org/lsst 

4 https:/ /www. darkenergysurvey.org/ 

5 http://jdcm.gsfc.nasa.gov/ 



signal can be measured by cross-correlating the CMB 
with tracers of LSS and has recently been detected us- 
ing WMAP data of CMB in combination with several 
LSS surveys at the the ~ 3 — 4cr confidence level pro- 
viding independent evidence for t he existenc e of th e DE 
(see iGiannantonio et all ()2008al ): iXia et all (|2009f ) and 
references therein). 

Cross-correlation then provides a powerful method to 
discriminate among different DE models and, in particu- 
lar, to detect a possible interaction between DE and DM 
other than investigate the clustering properties of DE on 
large scales. If present, DM-DE coupling changes both 
the scaling of the DM energy density and the growth rate 
of matter perturbations yielding a significant evolution of 
the metric potentials even in the matter era. 

The aim of this paper is to provide an exhaustive anal- 
ysis of the ISW effect in the contest of t he so called 
coupled DE cosmologies ([Amendola 1 120001 ) mainly fo- 
cusing on the effects of the DM-DE coupling on the 
cross-correlation between galaxy/quasar distributions 
and ISW effect. 

Such models can be motivated in the contest of scalar- 
tensor theories of gravity or describe the low energy limit 
of a more fundamental theory beyond the standard model 
of particle physics, e.g. string theory. 

The models which we aim to investigate differ from the 
standard ACDM in three different aspects: (i) DE is a 
self-interacting scalar field <j> rather than a cosmological 
constant A. We shall consider a class of self-interaction 
potentials V(<f>) admitting tracker solutions, (ii) A linear 
DM-DE coupling is allowed, (iii) We allow neutrinos to 
be massive. 

The effects of massive neutrinos in cosmology have 
been studied thoroughl y for m any years (for a review see 
iLesgourgues fc Pastor I (|2006l )). Cosmological observa- 
tions are mostly sensitive to the sum of neutrino masses, 
M u . Currently, the strongest upper limit on neutrino 
mass scale comes from cosmology. One of the effects of 
massive neutrinos is to induce a small decay of the grav- 
itational potentials during both matter and DE domi- 
nation so that, in principle, ISW effect would also pro- 
vide inform ation on their m a ss. Fur the rmore, as recently 
outlin ed in lLa Vacca et al 1 (2009a. b); iKristiansen et al I 
(2010) the effects that massive neutrinos have on the an- 
gular power spectrum of the CMB anisotropics, C;, and 
matter power spectrum, P{k) 1 are almost opposite to 
those of the DM-DE coupling, resulting in a strong de- 
generacy between the coupling strength /? and M v . A 
recent analysis by means of Monte Carlo Markov Chain 
method has shown that a cosmology with significant M„ 
and /? is statistically preferred to one with no coupling 
and almost massless neutrinos. Further, when priors on 
the neutrino mass from earth-based neutrino mass exper- 
iments (Heidelberg-Moscow neutrinoless double /3-decay, 
KATRIN tritium /3-decay0) are added to the analysis, a 
5 — 6cr detection of a DM-DE coupling is found. 

The plan of the paper is as follows: in Section 2 we 
describe our model while the ISW effect theory is re- 
viewed in Section 3 where we also discuss how the ISW 
signal depends on the main parameters of the model. In 
Section 4 we discuss galaxy bias and magnification bias. 

6 http: //www-ik. fzk.de/katrin/publications/documents/Design 
Report2004-12Jan2005.pdf 
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Comparison between theoretical prediction and observa- 
tion data is presented in Section 5 while a tomographic 
analysis is performed in Section 6. Section 7 is devoted 
to the conclusions. 

2. THE MODEL 

We assume a spatially flat Friedmann-Robertson- 
Walker (FRW) background with metric ds 2 = 
a 2 (rj) (—drj 2 + dx l dxi) (rj is the conformal time) filled 
with baryons, photons, neutrinos, DM and a component 
of DE which will be ascribed to a scalar field <j> self- 
interacting through a potential V(4>). In the following, 
the indexes 6, c, v and <f> will denote baryons, cold DM, 
massive neutrinos and DE. Photons and massless neutri- 
nos will be referred as radiation and denoted by r. 

In addition to self-interaction we also consider a possi- 
ble interaction between the scalar field and DM. Here we 
give only the equations for baryons, DM and DE being 
the equ ations for the other compo nents the usual ones 
(see, e.g lMa fc Bertschingerl (|1995l )b 

The Friedmann equation for the scale factor a, the con- 
tinuity equations for baryons and DM and the evolution 
equation for the scalar field read: 



n 2 = f G( Pb 



Pc + P<j>) a 
pb + SHpb = 



(1) 
(2) 
(3) 
(4) 



p c + WH p c = -C4>p c 

4> + 2H4> + a 2 V((t>) =Ca 2 p c 

where ' denotes the derivative with respect to rj, 
% = a/a, pi is the energy density of the component 
i = b, c, v, 4>, r and the constant C parametrizes the DM- 
DE coupling strength. 

Working in the conformal Newtonian gauge the metric 
of a perturbed flat FRW Universe takes the form: 

ds 2 = a 2 (r) [- (1 + 2$) dr 2 + (1 - 2*) dx l dxi\ (5) 

where <& plays the role of the Newtonian potential, \E' is 
the Newtonian spatial curvature, |$|, |\&| << 1 and: 



3 

2 



(6) 



where = Pi/p er , $i = <>Pi/pi, Qi, &i an d Wi = Pi/ Pi are 
the density parameter, density contrast, four-velocity di- 
vergence, shear and state parameter of the component i 
(pi and per — 3"H 2 /87ra 2 being the pressure of the com- 
ponent i and the critical density of the Universe). 

Linear perturbation equations for DM and baryons 
read: 

5 C + 9 C - 3* = -C6'<t> (8) 
8 c + (H-C<j>)8 c = k 2 (<S>-C5<j)) (9) 
6 b + e b - 3* = (10) 
6\+U9 b -c 2 s k 2 8 b = k 2 <$> + T phot _ b (11) 

where c 2 is the baryon sound speed, T p hot-b is the stan- 
dard term describing momentum exchange with photons 



due to Thomson scattering (see, e.g Ma & Bertschinger 
1995) and 5(f) is the perturbation to the scalar field which 
evolves according to: 



5<j) + 2U5<j) + {k 2 



a 2 V")5(f> - 4$0 + 2a 2 $U' = 

C(pA + 2p c $)a 2 (12) 

As we are interested in the cross-correlation between 
ISW effect and galaxy distributions, above equations can 
be simplified. At late time radiation and massive neutri- 
nos can be neglected so that no shear stresses are present 
and $ = ^. Furthermore, the cross-correlation signal 
comes from scales well within the horizon, ~ 100 — 200 
Mpc, so that the second term in (6) can be neglected. 

The above equations then reduce to the usual Poisson 
equation for the gravitation potential: 

3 V 2 

^ = ---^[n c s c + n b s b + n v 8 v } ) (13) 

a modified Jeans equation for DM and the usual one for 
baryons: 



5 C + {H- C<P)S C = -H 2 



l + ^(3 2 ) <>,A. <>;A - <>,A 



s b + ns b - -h 2 pcSc + n b 5 b + n„<y„] = o, (14) 

and a Poisson-like equation for the scalar field perturba- 
tion: 

T-L 2 

8(t> = -p- [n c s c + n b 5 b + Q, V 5 V ] (15) 

where we have defined the dimensionless coupling pa- 
rameter: 

P = V3/16vr m p C 

(m p — G -1 / 2 is the Planck mass). 

As clearly visible from the ab ove equations and widely 
discussed in lAmendola I (|2000D . coupling affects the dy- 
namics of DM particles. As a consequence baryons and 
DM develop a bias 6*, i.e. 6 b = b*6 dm . Notice that, 
this bias, which origin is to ascribe to the coupling, is 
something completely different from the galaxy bias due 
to hydrodynamical effects, discussed in the subsequent 
sections. 

It is also worth mentioning that, unlike the uncoupled 
case, in the presence of coupling, Universe goes through 
an evolutionary phase named </>-matter dominated era 
(0MDE) just after matter-radiation equivalence. In 
this period the scalar field <j) behaves as stiff matter 
{P4>/ Pij> = 1) having a non-negligible kinetic energy which 
dominates over the potential one. After this stage, the 
usual matter era follows before entering in the acceler- 
ated regime with a final De Sitter attractor. Notice also 
that, because of the 0MDE and the non-usual scaling of 
the DM energy density, i.e. p c cx a~ 3 e~ c ^, after equiva- 
lence the background expansion law will differ from the 
usual a (x rj 2 . 

2.1. Potential 

We shall consider the Ratra-Peebles (RP) 
dRatra fc Peebles I 119881) an d SUGRA self-interaction 
potentials ()Brax et al Il2000| ). reading 



A a+4 



RP (16) 



4 



+ 



1.2 
1 




p=0, m v =0 eV 

p-0, m v -0.15eV 

p-0, m v =0.46 eV 

P-0.075, m v -0 eV 

p=0.15, m v =0eV 

ACDM 


0.8 






0.6 


V 




0.4 


SUGRA 
k=0.027 h/Mpc 




0.2 







1e-05 0.0001 0.001 0.01 0.1 1 

a 

Fig. 1. — Evolution of the gravitational potentials as a function 
of the scale factor for different values of and m v for SUGRA 
model. For comparison the ACDM case is also displayed 



and 

A«+ 4 4tt4- 

V ^) = — e P SUGRA (17) 

respectively; they allow tracker solutions for any a > 0. 
For both potentials, once a and A are assigned, the 
present time DE density parameter flng is uniquely de- 
fined. 

Limits on these models without coupling between DE 
and CDM have been studied in La Vacca & Kristiansen 
2009. They find, that only quite small A = log A/GeV 
are allowed. In the SUGRA case in particular only 
A < —3.5 is allowed. Such small values are well below 
the range motivated by particle physics. Therefore the 
physical appeal of the SUGRA potential is spoiled. 

Let us however outline that, when the [3 degree of free- 
dom is opened, A values as large as 30 GeV become al- 
lowed, at the l-<7 level, while, at the 2-a level, no sig- 
nificant constraint on the energy scale A remains. Even 
for the RP potential, for which a limit A < —8.5 held, in 
the absence of coupling, values A <~ — 2 become allowed 
(La Vacca et al 2009). 

In the absence of DM-DE coupling, RP yields quite a 
slowly varying w^a) — p^j ' state parameter. On the 
contrary, SUGRA yields a fast varying w,p. Although 
coupling causes a behavior significantly different from 
the uncoupled case, one could again consider these po- 
tentials as examples of rapidly or slowly varying w^. 

The results shown in the next two sections are quali- 
tatively the same for RP and SUGRA models. We will 
show them only for SUGRA cosmologies while the re- 
sults for RP models will be shown only when comparing 
the theoretical predictions with observational data and 
dealing with redshift tomography. 

3. ISW EFFECT 

ISW effect arises when CMB photons from the last 
scattering surface pass through a time-dependent gravi- 
tational potential changing its energy so that additional 
temperature anisotropies are generated. The ISW tem- 
perature fluctuation, AT ISW , is given by: 

AT isw = T f e- T ($ + *)rfr? . (18) 



where T is the CMB temperature and e~ T is the visibility 
function of the photons. 

As outlined in previous section, we will deal with scales 
within the horizon and redshifts such that radiation and 
any anisotropic stress can be neglected ($ = *). In 
the matter era and in the absence of DM-DE coupling 
(P = 0), Poisson equation reads: 

3 H 2 

^ = -^-j^n m 5 m (19) 

(fl m and S m oc a cx nf are the total matter density pa- 
rameter and density contrast respectively) from which 
one can appreciate that the gravitational potential stays 
constant, <f> = 0, and no ISW effect arises. However, 
when DE starts to dominate the cosmic expansion, $ is 
no longer constant and ISW effect generates secondary 
anisotropies in the CMB. 

On the other hand, as explained in the previous sec- 
tion, DM-DE coupling affects both the background and 
density perturbation evolution, resulting in a variation 
of $ even during the matter domination. 

Figure 1 shows how [3 and m v affect the time evolu- 
tion of the sum $ + which time derivative forms the 
source of the ISW effect. Evolution of the gravitational 
potentials is obtained by a modified version of the public 
code CMBFAST integrating the fully relativistic equa- 
tions and taking into account the contributions of all of 
the components, i.e, photons, DM, baryons, neutrinos 
and DE. 

Notice how (3 and ra v affect $ + ^ in an opposite fash- 
ion. 

Performing a measurement of the ISW effect is, how- 
ever, a difficult task because of its small signal compared 
with that of primary anisotropies (~ 10 times larger). 
Furthermore, while on small scales the small differences 
in temperature tend to cancel out, the large scales, from 
which the most ISW effect contributes come from, are 
strongly affected by the cosmic variance. 

The problem can be overcome by cross-correlating the 
ISW anisotropies with some tracers of the matter density, 
e.g. astrophysical objects like galaxies. 

The observed galaxy density contrast in the direction 
rii is: 

S gal (n 1 )= Jb{z)^(z)5 m (n lt z)dz (20) 

where dN/dz is the normalized selection function of the 
galaxy survey and b(z) is the galaxy bias, which will be 
discussed in the next section, relating the galaxy den- 
sity contrast to the inhomogeneities in the mass distri- 
bution, 5 g a i = bS m . Since 6 m is related to the gravita- 
tional potential through the Poisson equation, the ob- 
served galaxy density will be correlated with the ISW 
anisotropies in the nearby direction n 2 : 

AT ISW (n 2 ) =2T J e- T ^^(n 2 ,z)dz (21) 

The 2-point angular cross-correlation function (CCF) 
and auto-correlation function (ACF) in the harmonic 
space are then defined as: 

C ISW -^ l (e) = (AT(n 1 )S gal (n 2 )) 
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Fig. 2. — Redshift evolution of $ + ^ (left) and its time derivative (right) for different values of A in uncoupled SUGRA with masslcss 
neutrinos 





Fig. 3. — ISW-matter cross-correlation power spectra (top) and functions (bottom). Their dependence on A is shown at z=0.3 (left) and 
z=3 (right). Models are the uncoupled SUGRA with massless neutrinos 
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C^-^(0) = (<W(Ai)V(n 2 )) 

= Y,-^r C r l - 9al Pi[cos(6)} (22) 



1=2 



where 6 = |ni — A2I, Pi's are the Legendre polynomials 
and the cross- and auto-correlation power spectra arc 
given by: 



a 



ISW-gal 



J d ^\k)ii sw {k)ir\k) 



where A 2 is the primordial power spectrum of scalar per- 
turbations and the integrands I ISW and I gal are: 

ll sw (k) = 2T J e- T ^ji[k X {z)]dz , (24) 
ir\k)= Jb(z)^(z)5 m (k,z)ji[k X (z)}dz (25) 

(here ji(x) are the spherical Bessel functions, and x 1S 
the comoving distance). 

In the following we use our modified CMBFAST code 
to calculate the theoretical CCF and ACF. In order to 
better understand the effects of A, (3 and ra„ on them, 
we start to compute the ISW-matter CCF and ACF, 
C ISW - m {6) and C m - m (<9), and power spectra, c\ sw - m 
and C ; m ~" 1 , the values of which are obtained similar to 
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Fig. 4. — Redshift evolution of $ + (left) and its time derivative (right). Models are coupled SUGRA with massless neutrinos 
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Fig. 5. — ISW-matter cross— correlation power spectra (top) and functions (bottom). Their dependence on A is shown at z=0.3 (left) and 
z=3 (right). Models are coupled SUGRA with massless neutrinos 



and (f2"3"| by replacing 6 ga i(hi) and If al {k) with 



8 m (hi) = J ^-(z)S m (h 1 ,z)dz 



ir(k) 



dN 

dz 



(z)S m (k,z)ji [kx(z)]dz 



(26) 
(27) 



where the only difference from 5 ga i and If al is the bias 
factor. Here, we model dN/dz as a narrow Gaussian 
centered at two different redshifts, z=0.3 and z=3. This 
will permit to obtain some information about the time 
evolution of the correlations. Realistic selection functions 
will be considered in the next sections. Cosmological 
parameters are assumed to be t he same as in the WM AP 
5 years best fit ACDM model (jKomatsu et al II2009D . 

3.1. Dependence on A, (3 and m„ 



Let us consider first the uncoupled case ((3 = 0) and no 
massive neutrinos. Figure [5] shows the redshift evolution 
of the gravitational potentials (left panel) and the source 
of ISW effect $ + (right panel) for different values of 
A. The dependence on A of the ISW-matter correlations 
c iSW-m and c isw-m^ is then displayed in Fig. [3] 
at two different redshifts, z = 0.3 and z = 3. When 
increasing A, both c{ sw ~ m and C ISW ~ m (6) show an 
opposite behavior at low and high redshifts. This reflects 
the behavior of $ + 4'. 

In the presence of coupling one can distinguish be- 
tween two different behaviors for small and large A's. 
This is shown in Figs. |4] and [5] for j3 = 0.1. In the 
first case, the evolution of the gravitational potentials 
and the cross-correlation signal are almost independent 
from A. It can be understood noticing that for small A, 
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Fig. 7. — ISW— matter cross-correlation power spectra (top) and functions (bottom). Their dependence on /3 is shown at z=0.3 (left) and 
z=3 (right). Models are the SUGRA models with massless neutrinos 



0MDE is very long and the usual tracker solution is (al- 
most) never reached. In this phase, coupling terms in the 
field equations dominate so that the tracker solution is 
almost independent from A. When increasing A, </>MDE 
becomes shorter and the behavior resemble that of the 
uncoupled case. The transition between the two above 
regimes occurs around A = 1 GeV. 

Dependence on /3 is shown in Figs |6] and [7] Again, 
the behavior of the cross-correlation reflects that of the 
ISW source. However, while coupling can have opposite 
effects on the cross-correlation at different redshifts, i.e. 
it can increase or decrease the signal, massive neutrinos 
always decrease the signal. It is shown in Fig. |8] which 



displays the behavior of C x 



ISW-m 



as a function of m„ at 



two different redshifts. Cross-correlation signal always 
decreases at the increasing of the neutrino mass. 



4. GALAXY BIAS AND MAGNIFICATION BIAS 

The galaxy bias b can, in general, evolve both in red- 
shift or as a function of the scale. However, on the large 
scales of interest for the ISW effect, the bias is usually as- 
sumed to be linear, spatially constant and only redshift 
dependent, i.e. S ga i = b(z)S m . This assumption is fully 
consistent with results from numerical simulations, red- 
shift surveys and semi-analytic ca lculation in the contes t 
of the so-called halo -model (see iBlanton et all (|1999ft ; 
IPercival et all ([2Q07f )). 

However, given a galaxy selection function dN/dz 
picked at certain redshift z, the bias can be approxi- 
mated with an appropriated constant. In this case it will 
be 

£jgal — gal j^2^j7n— m (28) 

Within the above approximation, given a particular sur- 
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Fig. 8. — Effect of massive neutrinos on ISW-matter cross cor- 
relation 



Fig. 9. — Comparison between cross-correlation spectra with and 
without magnification bias at z=1.5 and z=3.0 



vey and assumed a cosmological model, the bias is usu- 
ally estimated by fitting the theoretical matter-matter 
correlation function, C m ~ m , for the assumed cosmol- 
ogy, to the observed galaxy-galaxy correlation function, 

^jgal — gal 

Biases have been estimated for different surveys 

by several authors ass uming the WMAP b est fit 

ACD M cosmology Csee [Boughn fc Critteden J (12002L 
2004ft: iGiannan tonio ct al I (I2006D: iMyers etaTI (|2006t ); 
Rassat et al I (|2007ft : iBlake et al I (|2007ft ) Since we are 
considering cosmological models different from a ACDM, 
we need to appropriately rescale those biases to each of 
the our models. Note, however, how the estimation of b in 
(|28[) depends on the normalization of the power spectrum 
in C m ~ m (see (f2"3")l ). For a fixed normalization, taking 
into account (f25|) . biases will be rescaled according to: 

. *~tm—m . 

"model — °ACDM r m-m V /y ) 

^ ^model > 

where < > indicates the average on the angular scales 
9 of interest. 

However, (|29| should be generalized when magnifi- 
cation bias effect due to gravitational lensing is non- 
negligible. Gravitational lensing by intervening matter 
changes the observed galaxy number density 8 ga i , leading 
a correction term 5^ being added to the intrinsic galaxy 
fluctuation 5 ga i 

3gal = bgal + 5^ 

With this correction, the observed ACF becomes: 

^jgal — gal ^jgal — gal , ^^jgal — , ^j^l—^l 

= b 2 C m - m + 2bC m ~ fl + C"^ (30) 

where C x ~ v =< 5 x S y > and the rescaled bias will then 
be the solution of: 

°model <■ ^model > + ZD model < ^ model > "+" 

<C m 7del>~<C 9 A a CD 9 M>-0 (31) 

Auto- and cross-correlations corrected for magnifica- 
tion bias are ob tained considering in (25) the function 
(|Ho et al 112001 : 
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f{z) = b(z)^—6 m (k 7 z) + / dz'W(z, z'){a{z') - 1)^1 

(32) 

where a(z') is the slope of the number counts of galaxy 
number density as a function of the flux, A^(> F) oc 
F~ a . It depends on the choice of galaxy sample and is 
redshift dependent. The lensing window function (in a 
flat Universe) is: 

W(z,z') ^ kH(k,z)^^l X (z) 

Magnification bias increases with redshift and could be 
important when dealing with deep survey, e.g. quasars. 
This is shown in Fig. |9] which compares the effect of 
the magnification bias on the ISW-gal correlation at 
z = 1.5 and z = 3. A detailed analysis on how mag- 
nification bias affects ACF and CCF c an be found in 
lLoVerde et all (12001: iHui et al I (pUol : lLoVerde et ail 
<|2008f k lHui et all (|2008h .~ 

5. COMPARISON TO OBSERVATIONAL DATA 

Investigations of CMB-LSS correlations were made 
in a recent series of works which rely on WMAP 
data and a variety of LSS probes (iNolta et al l 
2001 lAfehordi et al I [20(51 iCabre' et al I 120061 1 20071: 
Rassat et al 20071: IRaccanelli et al ll2008HHo et al 1 12008: 
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Fig. 10. — Redshift distribution dN/dz of all catalogues 
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Fig. 11. — The observed CCF of six different galaxy catalogues. The curves are the theoretical predictions for the best fit ACDM, SUGRA 
and RP cosmologies (see text) 



Parameter 


RP 


SUGRA 


10 2 uj b 


2.260 ± 0.061 


2.260 ± 0.065 


Ulc 


0.1039 ± 0.0062 


0.1042 ± 0.0084 


T 


0.087 ± 0.016 


0.088 ± 0.017 


M„ (eV) (95% C.L.) 


< 1.13 


< 1.17 


(95% C.L.) 


<0.17 


<0.18 


log 10 (A/GeV)(95%C.L.) 


< -4.2 


< 6.3 


n a 


0.969 ± 0.015 


0.970 ± 0.018 


ln(10 10 A s ) 


3.055 ± 0.040 


3.057 ± 0.041 


H (km/s/Mpc) 


71.8 ± 2.5 


71.9 ± 2.7 



TABLE 1 

Best fit values and 1 — <t error bars for RP and SUGRA 
models. Only upper limits on M„, /3 and A are shown. 





ACDM 


RP 


SUGRA 


2MASS 


1.40 


1.46 


1.47 


SDSS gal 


1.00 


1.03 


1.04 


LRG 


1.80 


1.83 


1.84 


NVSS 


1.50 


1.53 


1.54 


HEAO 


1.06 


1.09 


1.09 


QSO 


2.30 


2.33 


2.33 



TABLE 2 

Galaxies biases for different catalogues and models. 
Biases are calculated a ccording to i|29[i using for & a(7om 

THE VALUES GIVEN IN Gl ANN ANTONIO ET AL ( 2008aD. 
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Fig. 12. — Splitting schemes (i) (left panel), (ii) (middle panel) and (iii) (right panel) described in the text. 



iGiannantonio et al |[2008at Kia et al 112009ft . There is an 
overall agreement among different groups in finding an 
evidence for a positive ISW effect at the ~ 3 — 4c con- 
fidence level. It was also found a substantial agreement 
between the observed cross-correlations and the expected 
signal arising from the ISW effect in the WMAP best fit 
ACDM cosmology. Differ ent DE mod 
els we r e also considered inlOlivares et al I (|2008ft ; iSchafer I 
(12008ft ; IGiannantonio et al I (l2008bl) . 

In a recent work ([Giannantonio et al I I2008aft a com- 
bined analysis of the ISW effect was performed by cross- 
correlating CMB map provided by the WMAP collabo- 
ration with all the relevant large scale data sets and mod- 
cling their covariance properties with different methods. 

In this section, we compare our theoretical predictions 
based on the m odels described abov e with t he measure- 
ments made in IGiannantonio et al I (|2008aft by consid- 
ering six different galaxy catalogues: the optical Sloan 
Digital Sky Survey (SDSS), the infrared 2 Micron All- 
Sky Survey (2MASS), the X-ray catalogue from the High 
Energy Astrophysical Observatory (HEAO) and the ra- 
dio galaxy catalogue from the NRAO VLA Sky Survey 
(NVSS). In addition, given the high quality of the SDSS 
data, some further subsamples was extracted from it, 
consisting of Luminous Red Galaxies (LRG) and quasars 
(QSO). 

As shown in Fig. [10] their redshift distributions, 
dN/dz, span a redshift range < z < 2.5. In order of 
increasing mean redshift of the sample we have: 2MASS, 
SDSS galaxies, LRG, NVS S, HEAO and QSO (for de - 
tails on the cat alogues see [Giannant onio et al I (|2008al) ; 
IHo et al I ((2008ft . 

Constraints from CMB, SNIa and LSS on cosmologies 
wi th coupling and massive neutrinos have been obtai ned 
in ILa Vacca etal! (|2009bf ): IKristiansen et al I (|2010ft by 
means of Monte Carlo Markov Chain technique. Best fit 
parameters from their analysis will be used in this and 
in the next section when dealing with redshift tomogra- 
phy. Parameters are summarized in Table [T] where cj&. c 
are the physical baryon and cold dark matter density 
parameters, ojf,. c = flb iC h 2 , where h is the dimensionless 
Hubble parameter, r is the optical depth to reionization, 
n s is the scalar spectral index, A s denotes the amplitude 
of the scalar fluctuations at a scale of k = 0.05Mpc _1 , 
M v = Em,,, assuming 3 equal neutrino masses m v , A 
denotes the energy scale in DE potentials, while (3 is the 
coupling parameter between DM and DE. In the follow- 
ing we will use A = lO^GeU and A = lGeV for RP 
and SUGRA respectively which correspond to the ^For 



the three splitting schemes described above 1-a limits. 
Best fit parameters for coupling and neutrino mass are 
approximately the same for both model, i.e. /3 ~ 0.1 and 

n v = M u /h 2 m.uev ~ o.oi. 

For each catalogue and model, biases are shown in Ta- 
ble 1 and calculate d according to (l29l using for &acdm 
the values given in IGiannantonio et al I (|2008aft . Given 
the low mean redshifts of the catalogues we neglect the 
magnification bias effect which amounts to a few percent 
only in the case of quasars. It will be, however, con- 
sidered in the next Section when dealing with redshift 
tomography and higher z. We will also discuss how good 
to approximate b with a constat is. 

For each catalogue, we then determine the expected 
CCFs for our models. 

Comparison with observational data is shown in Fig. 
ITTI for SUGRA and RP models. The predictions for the 
ACDM is also displayed. Note that, because of know con- 
tamin ation from Sunyaev- Zeldovich effect in the 2MASS 
data (jAfshordi et alll2004ft , the four smallest angle bins 
should be disregarded. While it is not possible to dis- 
tinguish among the models at low redshifts, discrepan- 
cies between coupled models and ACDM increase with 
z even though RP and SUGRA models remain indistin- 
guishable. In spite of this, however, current data alone 
seem not able to discriminate between coupled models 
and ACDM. 



6. REDSHIFT TOMOGRAPHY 

As already outlined, unlike uncoupled DE models with 
massless neutrinos, both coupling and massive neutrinos 
causes the gravitational potentials to evolve even in the 
matter dominated epoch. Therefore, a detection of a 
non-vanishing ISW effect signal at such high redshifts 
would rule out a vast class of DE models indicating a 
possible interaction in the dark sector. Upcoming galaxy 
surveys will cover a large redshift range. One goal will 
be to use the photometric redshifts of the galaxies to 
split the survey into multiple redshift bins allowing for 
tomographic analysis. 

Following the procedure of IHu fc Scranton I (|2004ft . 
given a galaxy distribution, n{z) = dN/dz, the galax- 
ies can be divided into photometric bins, labelled with 
index i: 
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Fig. 13. — Relative error on snr in the i-th tomographic bin for the cases in Tableland the three different splitting schemes described 
in the text. 




Fig. 14. — Same as Fig. 1131 but for the cumulative snr. 
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Fig. 15. — Relative error on snr t (top panel) and 
(bottom panel) for the cases in Tabl^S] We assume a shot noise as 
for the SDSS DR6 survey and compare with the case of negligible 
shot noise. 



Assuming a distribution n(z) of the standard form: 



n(z) 



p/m+l\ m 



-exp 



(33) 



and that the photometric redshift errors are Gaussian 
distributed with an rms fluctuation er(z), the resulting 
photometric redshift distributions are given by: 



ni(z) = \n{z) 



erfc 



Zi-i - z\ ( Zi — z 

— erjc 



V2a(z) 



\V2a(2 



We now propose to study how a tomographic analysis 
is affected when considering different splitting schemes. 
This will permit to single out the optimal splitting 
choices which guarantee a signal to noise ratio, snr, high 
enough to distinguish among different DE models. 

As we are interested in high z, we model the overall 
distribution n{z) according to (|33p with m = 2, f3 = 2.2 
and zq = 1.62. These values provide a good fit of the 
quasar distributi on from SPSS DR6 considered in the 
previous section (jXia et al 1120091 ) . We assume the shape 
of such a distribution to be approximately the same as 
for that expected from future surveys. 

Further, we assume a(z) = 0.03(l+z) as expected from 
future experiments, and consider three different splitting 
schemes, each with 5 bins in the redshift range from z = 
0.75 to z = 4: 

(i) bins equally spaced in z, with Az = Zi—Zi-i = 0.65; 

(ii) same number of galaxies, An, in each bin; 

(iii) bin sizes increasing proportionally to the photo- 
metric error, Az cx cr(z). 

The three splitting schemes are shown in Fig. [T2l 
The thick line is the overall quasar distribution while 
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Fig. 16. — The splitting scheme (i) (see text) is shown in the top panel on the left. Other panels show the cross-correlation signal in the 
five bins considered. Lower frames of each panel display the ratio between the SUGRA and RP spectra and the ACDM spectrum 



the other curves are the true (spectroscopic) distribu- 
tions that correspond to the divisions (vertical lines) in 
photo-z space. 

We also take into account magnification bias effect 
which can be importan t at higher z. In a recent anal- 
ysis bv IHo et all (|2008f ). the slope of the quasar counts, 
a, entering in (|32j) . was found to be redshift dependent. 
They found a — 0.82 in the photometric redshift range 
0.65 < Zp ho to < 1-45 and a — 0.9 for 1.45 < z photo < 2. 
For simplicity, w e assume a constant slope a = 0.9 as in 
IXia et al I ff2009h . 

6.1. Dependence on galaxy bias evolution and 
magnification bias 

Before comparing the different models we discuss the 
effects of ignoring magnification bias and quasar bias evo- 
lution. This is done for the ACDM cosmology. Results, 



however, are valid for RP and SUGRA as well. 

Notice that implications of galaxy bias evolution on 
ISW measurements and parameter estimation has been 
considered in a pioneer work of Schaefer et al 2009. 

For the quasar bias evolu tion in ACDM, we use the 
empirical formula derived bv iCroom et al I ((2004): 

b(z) = 0.53 + 0.289(1 + z) 2 (34) 

which p rovides a good fit of the recent observational find- 
ings bv lXia et all (120091 ). 

We first consider the overall quasar distribution and 
calculate the expected signal-to-noise ratio, snr, of the 
cross-correlation in the following cases: constant quasar 
bias b — b{z) (where z is the mean redshift of the survey) 
and b = b(z) as given by (|34|) . For each of them, the snr 
is obtain by neglecting or considering the magnification 
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Fig. 17. — The expected ISW-LSS cross correlation signal to noise ratio, snr, for the different splitting schemes considered in the text 
as a function of the mean redshift of the bins for the best fit ACDM, RP and SUGRA cosmologies. 
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Fig. 18. — Same as previous Figure but for the cumulative snr. Horizontal lines indicate the values of snr obtained by using the overall 
distribution. 
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Fig. 19.— RP and SUGRA snr's of Fig. [17] (top panel) and Fig. 
1181 (bottom panel) normalized by the ACDM ones, snr obtained 
by using the full distribution are also shown. 



bias correction. Results are summarized in table [3] 

F or Gaussian file ds, the expected snr is given by (see 
e.g. ICooravl (pOOl O: 

U.« 

(snr) 2 = f sky J2 ( 2l + 1) x 

Imin 



quasar bias 


magnification bias 


snr 


b = b(z) 


no 


4.56 


b = b(z) 


no 


4.31 


b = b(z) 


yes 


4.61 


b = b(z) 


yes 


4.38 



TABLE 3 

Effects of ignoring galaxy bias evolution and 
magnification bias on the snr for acdm. z is the mean 
redshift of quasar distribution used (see text). 



quasar bias magnification bias 



a) 


b t 


= b{z,) 


no 


b) 


bi 


= b(z) 


no 


c) 


b, 


= b(zi) 


yes 


d) 


b t 


= b(z) 


yes 


e) 


h 


= b(z) 


yes 



TABLE 4 

Cases considered in order to discuss the effect of 
ignoring galaxy bias evolution and magnification bias on 
snr when dealing with tomography. zi and z are the mean 
redshift of the i-th bin and of the overall distribution. 
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Fig. 20. — Left panel: expected snr values for the best fit RP model and splitting scheme (iii) with 2, 3,5 10 bins, snr obtained by using 
the full distribution is also shown. Right panel: snr values of the left panel normalized by those expected in the ACDM case. 



^plSW-galy 
f c ISW-gaiy + ^isw + C Nisw^ (pgal + 

(35) 

where C/ V/Sw = Cf + Cf et is the noise contribution to 
the ISW (Cf and Cf et being the total anisotropy con- 
tribution and any detector noise contribution, negligible 

on the scales we are interested in) while cf aal — 1/N, is 
the shot noise associated with the galaxy/quasar catalog 
(N is the surface density of galaxies/quasars per stera- 
dians). f s k y is the fraction of sky common to CMB and 
galaxy/quasar survey maps. 

Note that (f3"5|) is strictly true only in the case of Gaus- 
sian fields and full-sky coverage, f s k y = 1. In this sec- 
tion, we assume that to be the case other than a neg- 
ligible C^ gal . These assumptions will be relaxed in the 
following. Anyway, for partial sky coverage, one can, in 
first approximation, multiply for f sky the values of snr 
here presented. 

The value of the lowest multipole, l m im can be approx- 
imately set to l m in = 7T /2f s k y in order to account for the 
loss of low multipole modes. Anyway, although a signifi- 
cant part of the ISW signal comes from lower multipoles, 
we set l m in = 10 in order to avoid effects of ga uge correc- 
tion o n very large scales recently discussed by lYoo et al~1 
(2009) even if this implies a reduction of snr. The maxi- 
mum multipole, l m ax, is set to l m ax = 1000. However, we 
will show later that the contribution to snr from I > 400 



multipoles is negligible. 

From Table 03 it follows that cross-correlation mea- 
surements are more affected by errors when ignoring the 
quasar bias evolution rather than the magnification bias 
correction. In fact, in this last case, an error of ~ 1.6% 
on snr is obtained while in this last case while the error 
raises to ~ 7% if the quasar bias is approximated by a 
constat value, b — b(z) (z being the mean redshift of the 
quasar distribution) and it is of ~ 4% if both quasar bias 
evolution and magnification bias are neglected. 

We now turn to tomography. According to ([35|, the 
signal to noise ratios, snri, in the i-th tomographic bin 
are calculated for each of the cases listed in Table SJ 
Fig. [TS] shows the errors Asnr^ relative to the case e) for 
each bin of the three splitting schemes described above. 
Unlike the overall quasar distribution, in each bin, dis- 
tributions are very narrow and the quasar bias can be 
approximated by a constant, 6, = b(zi) (zi being the 
mean redshift of the i-th bin), leading to only a minor 
error < 2% on snri (case c)). On the other hand, ig- 
noring magnification bias correction might be critical at 
high z (Asnri/snri ~ 7.5%, case b)) and if a constant 
bias, b — b(z), is used for all bins, the error can reach 
~ 10% at higher redshifts. 

However, since photometric redshift errors cause the 
bins to overlap and magnification of the galaxies at Zi 
probes structures at z < z i: cross-correlation measure- 
ments at high z are quite correlated with those at low 
redshifts. Taking into account such correlations, the net 
accumulated snr for measurements from all the bins up 
to z max is given by: 



(snr(w)) 2 = E E Gl sw -^ [Covf\ C ; SW ~^ (36) 

Zi,Zj<Z max l m i„ 



(jlSW-gali QlSW-galj {qISW i qN isw \ ( fjgah ggalj ^ ^AgoiA 

[Govt]. , = ^ — (37) 

[f aky (2l + l)}- 1 



I 

The relative error on snr(z max ) is shown in Fig. [TJ] for the same cases as in Fig. 1131 In all cases, despite the 
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Fig. 21. — Cumulative snr from all the bins of the splitting scheme (iii) for ACDM, RP and SUGRA models. Left panel: contribution 
to snr from each multipole I. Right panel: contribution to snr up to I = l m ax- 



ACDM 





















^^^\^^^ 


:T" 


~~ 3 


IfcV ^^^^ 







50 100 
N / N SDSS 





Fig. 22. — Cumulative snr/Al contour levels in the plane f a ^ y — N /N$DSS- Here NgDSS is the number density of quasars from the 
SDSS DR6 survey. 



significant differences in snri in high z bins, errors are 
always < 2%. This is understood because most of the cu- 
mulative snr comes from low redshifts. However, in the 
case e) , errors considerably increase if the shot noise asso- 
ciated to the catalog is non-negligible. Fig. [TBI compares 
the relative errors on snri (top panel) and snr(z ma x) 
(bottom panel) in the case of negligible shot noise and 
assuming a SDSS DR6-like survey with a number den- 
sity of quasars, N, of ~ 120/deg 2 . In this last case, the 
error on the cumulative snr ranges from 20%, to 40%, 
depending on the number of bins considered while the 
error on snri can also reach the 50% in bins at higher 
redshift. 

6.2. Model comparison 

In the previous section we have shown that when deal- 
ing with narrow tomographic bins, we can approximate 
the galaxy bias in each bin with a constant committing 
an error of no more than a few percent. 

Then, after calculating the cross-correlation power 
spectra for the best fit ACDM, SUGRA and RP models 
we apply (f5T|) to each bin of the three splitting schemes. 
The rescaled biases b{zi)np and b{z,i)suGRA are then fit- 
ted with an expression similar to (|34[) obtaining: 

b(z) = 0.54 + 0.291(1 + zf 

for RP and 

b(z) = 0.54 + 0.287(1 + zf 



for SUGRA. These expressions, which will be used in the 
rest of the paper, are valid up to z = 4 leading a galaxy 
bias evolution only slightly different from that of ACDM. 

In Fig. 1161 we show the cross-correlation power spectra 
for the splitting (i). Lower frames of each panel display 
the ratio between the SUGRA and RP spectra and the 
ACDM spectrum. Mean redshifts of the true bin distri- 
butions are also indicated. A similar redshift evolution 
could be obtained by considering (ii) or (iii). As clearly 
visible from the figure, a better discrimination among the 
models is expected at higher redshifts. 

The expected snr,; and snr(z max ) for the different 
splitting schemes are shown in Figs. [17] and [18] respec- 
tively, for the cosmologies considered. Despite the same 
qualitative behavior for all the models, higher snr are ex- 
pected in the SUGRA case; RP being in the between of 
SUGRA and ACDM. Horizontal lines in Fig. IT8l indicate 
the snr obtained by using the full distribution. 

For each splitting scheme, in Figs. [HI we plot the ratio 
between the RP and SUGRA snri and snr(z max ) values 
and those expected in the ACDM case. A first thing to 
note is the overlapping, for both RP and SUGRA, of the 
three curves corresponding to the different splittings indi- 
cating that the three schemes considered perform equally 
well in discriminating among the models. A better dis- 
crimination is however achieved looking at high redshifts. 

In Fig. 1201 we investigate how snri changes when in- 
creasing the number of the bins (left panel) and whether 
a greater number of bins could permit to better discrim- 
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inate between models (right panel). Results are shown 
for the splitting scheme (iii) with 2,3,5,10 bins in the RP 
model. Similar results are obtained in the other cases. 
snr in the case of no splitting is also shown. More bins, in 
principle, would permit to have a more detailed descrip- 
tion of the redshift evolution of ISW effect. However, 
as clearly visible in the left panel, the snr decreases at 
the increasing of the bin number. In the right panel, 
RP model is compared to ACDM. Even though, at high 
rcdshifts, tomography permits to distinguish among the 
models better than using the full distribution, the figure 
shows that increasing the number of the bins from 2 to 
10 would permit only a minor improvement in discrimi- 
nating between the models. 

6.3. Sky coverage and shot noise 

Up to now, we have considered the ideal case of Gaus- 
sian fields, full sky coverage and negligible shot noise. 
However, after cutting out our galaxy from the analy- 
sis, future CMB and galaxy maps are expected to cover, 
at best, a sky fraction f s k y — 0.7 — 0.8. In this case 
different multipoles are no longer independent and (|35[) , 
(|36[) and (|37|) only provide approximated estimations and 
a more rigorous analysis taking int o account for the ef - 
fective survey geometry is needed (|Cabre' et al I (|2007h . 
iHivon et all (|2002t l. iXia et al I d2011j)). It has been how- 
ever shown that, under the above approximations, a bet- 
ter estimation can be obtained by binning the power 
spectra data in bins of appropriate size making the 
bins independent. In this case, (|35|) and (l36l) are in- 
creased by a multiplicative factor of AL iCabre' et al I 
(|2007h found that Al = 20, 16, 8, 1 works well for f sky = 
0.1,0.2,0.4,0.8. 

In Figs. |2"T1 and |2"21 we show some results for the cumu- 
lative snr from all the bins of the splitting (iii) . Results 
are however the same for the other schemes. 

The left panel of Fig. [21] shows the contribution to the 
cumulative snr from each multipole while the cumula- 
tive snr up to I — Imax is displayed in the right panel. 
The dependence on f s k y and Al has been removed. As 
clearly visible, most of the cross-correlation signal comes 
from lower multipoles and contributions from I > 400 are 
negligible. 

In Fig. 1221 cumulative snr contour levels are plotted in 
the plane f s k y — N/Nsdss where Nsdss is the quasar 
number density for a SDSS DR6-like survey. Given that 
such a survey cover ~ 20% of the sky, future experiments 
covering a sky fraction f s k y = 0.8 will increase the cumu- 
lative snr of a factor ~ 3 if the shot noise is reduced by 
1/10 and a factor ^3.5 — 4 in the case = IOONsdss- 
No significant improvement is obtained by further re- 
ducing the shot noise. For f s k y = 0.8 and negligible shot 
noise, the increasing of snr in the i-th tomographic bin 
can range from a factor of 4 (low z bins) up to 10 (high 
z bins). This is shown in Fig. [23] 

7. CONCLUSIONS 

In this work we have investigated ISW-LSS cross- 
correlation in coupled Dark Energy models with massive 
neutrinos. The presence of a coupling between DM and 
DE as well as massive neutrinos change both the back- 
ground and matter perturbation evolutions yielding, un- 
like the ACDM case, a time-variation of the gravitational 



potentials even during the matter domination. A signif- 
icant ISW signal is thus expected also at high redshifts. 

Firstly, we have investigated the dependence on the 
energetic scale, A, of the DE potential, the coupling 
strength j3 and the neutrino mass m v . We considered 
first the uncoupled case (/3 = 0) and massless neutrinos. 
We found that, when increasing A, both (jl sw ~ m and 
C ISW ~ m {6) show an opposite behavior at low and high 
redshifts. This in fact reflects the behavior of $ + \P. 
In the presence of coupling one can distinguish between 
two different behaviors for small and large A's. In the 
first case, the evolution of the gravitational potentials 
and the cross-correlation signal are almost independent 
from A. It can be understood noticing that for small 
A, coupling terms in the DE field equations dominate so 
that its solution is almost independent from A. When 
increasing A, the behavior resemble that of the uncou- 
pled case. Dependence on j3 was also investigated and, 
again, the behavior of the cross-correlation at different 
rcdshifts reflects that of the ISW source. However, while 
coupling can affect c\ sw - m (and C ISW ~ m {9)) in an op- 
posite fashion at high and low redshifts, massive neutri- 
nos always decrease the cross-correlation signal. 

Secondly, we have provided a simple expression, eq. 
(|29p . which permits to appropriately rescale the galaxy 
bias when comparing different cosmologies once the bias 
of a particular model, e.g. ACDM , is known and the nor- 
malization of the power spectrum in each model is fixed. 
We also give, a generalized version of (|29|) to the case 
when the magnification bias effect due to gravitational 
lensing is non-negligible (see eq. ([3T|) '). 

Then, we compare the theoretical prediction on the 
cross-correlation function for our models with the obser- 
vat ional data obtained for six d ifferent galaxy catalogues 
bv iGiannantonio et al I (|2008af ). We found that, while 
it is not possible to distinguish among the models at 
low redshifts, discrepancies between coupled models and 
ACDM increase with z even though RP and SUGRA 
models remain indistinguishable. In spite of this, how- 
ever, current data alone seem not able to discriminate 
between coupled models and ACDM. 

Finally, we studied the redshift tomography. Upcom- 
ing galaxy surveys will cover a large redshift range also 
providing photometric redshifts of the galaxies with high 
accuracy. This will permit to split a survey into multi- 
ple photometric redshift bins allowing for tomographic 
analysis. Here, we were interested to study how a to- 
mographic analysis of the ISW-LSS cross-correlation is 
affected when considering different splitting schemes and 
assuming photometric redshift errors as expected from 
future experiments. As we were interested in high red- 
shifts, where our models, unlike the ACDM case, are 
expected to provide a significant ISW effect signal, ISW 
effect was cross-correlated with quasars. The quasar dis- 
tribution was thus split in tomographic bins according to 
three different schemes: (i) bins equally spaced in z; (ii) 
same number of galaxies, in each bin; (iii) bin sizes in- 
creasing proportionally to the photometric error. Cross- 
correlation were then calculated in each bin. 

Our tomographic study was based on a signal-to- 
noise analysis. We started our discussion investigat- 
ing the effect on snr of ignoring the quasar bias evo- 
lution and magnification bias correction for an ideal sur- 
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Fig. 23. — Comparison between snri expected from a survey with 
fsky = 0.8 and negligible shot noise and the present one. The three 
splitting schemes are considered. 
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vey. We found that, if the overall quasar distribution 
is used, cross-correlation measurements are more af- 
fected by errors when ignoring the quasar bias evolution 
(Asnr/snr ~ 7%) rather than the magnification bias 
correction (Asnr/snr ~ 1.6%) 

However, when dealing with tomography the error on 
snri {i indicating the i-th bin) never overcome <~ 2.5% 
if the quasar bias in each bin is approximated with an 
appropriated constant, but it can reach ~ 7.5%, at high 
rcdshifts, when magnification bias is ignored. Errors on 
the cumulative snr, however, always stay below the 2%. 
On the other hand errors can increase up to 50% if the 
shot noise associated to the quasar survey is set to the 
current values. 

We then used the tomographic analysis in order to 
compare different cosmologies. We found that the above 
splitting schemes, perform equally well in discriminating 
among the models. A better discrimination is however 
achieved looking at high redshifts. 

We also investigated how the expected signal to noise 
ratio, snr, of the cross-correlation changes when increas- 
ing the number of the bins and whether a greater number 
of bins could permit to better discriminate between mod- 
els. Even though more bins would allow to have more 
information on the redshift evolution of the ISW effect, 
the snr decreases at the increasing of the bin number. 
As a consequence, although tomography, at high red- 
shifts, would permit to distinguish among the models 
better than using the full distribution, when comparing 
our models to ACDM it was shown that increasing the 
number of the bins from 2 to 10 would permit only a 
minor improvement in the discrimination. 

Finally, we showed that future wide field surveys 
{fsky ~ 0.8) can increase the cumulative snr of the cross- 
correlation of a factor <~ 3 (3.5 — 4) if the current shot 
noise is reduced by 1/10 (1/100) while the snr of the sin- 
gle bins can increase up to a factor 10 at high redshift. 

Our snr analysis suggest a discrimination power of 
future ISW-LSS cross-correlation measurements able to 
distinguish among different cosmologies. However, in or- 
der to assess the discrimination, more rigorous analysis in 
terms of Fisher Matrix and Monte Carlo Markov Chain 
are needed. They are currently under investigation and 
left for future works. 
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